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In the present work the power spectrum of a particular class of tachyon fields is compared with the 
one corresponding to a cosmological constant model. This is done for different barotropic indexes 70 
and the background space time is assumed to be of the spatially flat Friedmann-Robertson-Walker 
type. The differential equation describing the perturbations is solved numerically and the power 
spectrum at the scale factor value a = 1 is plotted for each case. The result is that the power 
spectrum of the standard tachyon field differs in many magnitude orders from the ACDM. However, 
the one with 70 = 1.91, which corresponds to a complementary tachyon field, coincides fairly well 
with the concordance model. Therefore, we conclude that the perturbed solutions constitute an 
effective method to distinguish between the different 70 values for the tachionization ACDM model 
and the fiducial model. The Statefinder parameters {r, s}, measuring the deviations of the analysed 
model from the concordance model, are also explicitly calculated. Our result suggest that, depending 
on the value of 70 , these models can explain the observed expansion history or the perturbation power 
spectrum of the universe, but they may have problems in describing both features simultaneously. 

PACS numbers: 98.80.-k; 98.80.Cq. 
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I. INTRODUCTION 

The vast majority of cosmological tests are in agree¬ 
ment with the simplest type of dark energy density, the 
cosmological constant, indicated with A, so the standard 
cosmology is usually denominated as ACDM model or 
the concordance model. Despite its striking observational 
success, the present value of the cosmological constant is 
largely unappealing from the theoretical point of view. 
The ACDM suffers from a variety of problems that range 
from quantum gravity considerations to fine tuning and 
coincidence problems mm- 

The standard concordance model is well studied. The 
formation of dark matter structures for this model is rel¬ 
atively well understood, thanks to both the linear per¬ 
turbation theory and numerical cosmological simulations 

i-lll- 

Cosmological simulations are one of the most powerful 
tools for the research of the non-linear evolution of large- 
scale structures of the Universe. The study of the power 
spectrum gives a wealth of information about the matter 
density held, and they are rehected in several measure¬ 
ments such as cosmological weak lensing. 

In the recent years some attention to tachyon driven 
cosmologies was paid [5]. Some reasons are that such 
cosmologies may play an important role in inhationary 
models as well as in the present accelerated expansion, 
simulating the ehect of the dark energy, E-m, depend- 
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ing upon the form of the tachyon potential m-m, m- 
[E]. These scenarios are inspired in some string cosmol¬ 
ogy models pp] . 

In references m and [22] there are considered three 
classes of tachyons, which are classihed in terms of the 
value of their barotropic index 70 (see below). On the 
same way, in [ 16 j it is shown that, in the limiting case of 
small A the standard and complementary tachyon helds 
dehned in that reference tend to the ACDM model. 

The tachyonization cosmological models described 
above constitute a large class of scenarios containing the 
ACMD as a particular limit. It is natural to search for 
observable ehects allowing to distinguish them from the 
ACMD. One of such effects is the power spectrum of the 
gravitational potential of the primordial fluctuations for 
tachyon cosmologies. These perturbations are responsi¬ 
ble for the large structure scale of the universe and, due 
to the non linear nature of the Einstein equations, an 
appreciable deviation between the two models may be 
observed. 

An important cosmological diagnostic pair {r, s} called 
Statefinder, was introduced by [23]. The Statefinder is 
a geometrical diagnostic and allows us to characterize 
the properties of dark energy in a model independent 
manner. The Statefinder is dimensionless quantity con¬ 
structed in terms of the scale factor of the Universe and 
its time derivatives only. This cosmological diagnostic 
will be applied below to the model under consideration. 

The aim of the present work is to compare the power 
spectrum between the concordance and the taquion mod¬ 
els mentioned above. The organization goes as follows. In 
section 2 the main aspects of the power spectrum of the 
primordial fluctuations in Friedmann-Robertson-Walker 
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(FRW) cosmologies are briefly discussed, and the differ¬ 
ential equation that describes the perturbations is de¬ 
rived. In section 3 the three classes of tachyons in terms 
of their barotropic index is reviewed, together with their 
defining equations. In section 4 the power spectrum of 
primordial fluctuations for these models is characterized 
by numerical methods and it is compared to the ACMD 
model. Section 5 contains the discussion of the obtained 
results. In section 6 , we applied the stateflnder diagnostic 
pair (dubbed stateflnder) to these models. 


II. COSMOLOGICAL SCALAR 
PERTURBATIONS 

The small scalar perturbations in a spatially flat 
Friedmann-Robertson-Walker scenario in the conformal 
Newtonian gauge are represented by the line element [24] 

ds^ = a{ri) [(1 -I- 2^)dr]^ ~ (1 ~ 2'i’)Sijdx^dx^] , (I) 

where a{T]) is the scale factor as a function of the con¬ 
formal time Tj. The last is defined by dr] = dt/a, up to 
an additive constant. In this gauge, the variables 4) and 
'k are gauge invariant. In the situations for which the 
spatial part of the energy-momentum tensor is diagonal, 
it follows that $ = 4/ = 0, with cj) the Newtonian poten¬ 
tial. So, there remains only one free metric perturbation 
variable (j). As can be seen from Q, the gauge invariant 
quantities $ and 4/ can be interpreted as the amplitudes 
of the metric perturbations in the conformal-Newtonian 
coordinate system. 

In the following, a perfect fluid with an equation of 
state which, for each fluid species, takes a barotropic 
form p = (jQ — l)p will be considered. The anisotropic 
stress leads to nondiagonal space-space components of 
the energy-momentum tensor and vanishes for a perfect 
fluid. In this situation </> = i/f in Eq. Q, as remarked 
above. The general from of the classical equations which 
describe the evolution of small perturbations in a hydro- 
dynamical universe [25] . can be written in the following 
form 


solution of the form exp(—ik.x), the equation (j^ can 
be cast in the form 


da^ 


+ a 



where H = a{t)/a{t) is the Hubble parameter and the 
dot denotes the derivative with respect to time t. 


III. TACHYON FIELD EQUATIONS AND 
ACDM 

The tachyon held of m,m an scalar field (j) described 
by the addition of the following lagrangian 

L = u(</.)yr^a;:p^, (s) 

to the Einstein one with non zero cosmological constant 
A. Here V (</>) is the self interaction potential [TB] for the 
tachyon. The background energy density and pressure of 
the tachyon condensate, for a flat FRW cosmology, are 
given by 

P<p = / ^ > Pet, = -V ( 6 ) 

respectively. The Friedmann and the standard conserva¬ 
tion equations become 

= K + —Z=, (7) 

VI 

-I-3i70(I — -I- 

respectively, where H = d(t)/a(t) is the Hubble parame¬ 
ter and a{t) is the cosmic scale factor and the dot denotes 
derivative with respect to the time variable t. The equa¬ 
tion of state for each fluid specie takes a barotropic form 
P = (70 ~ 1)P! where 7 is the barotropic index. 


u” — = H. ( 2 ) 

In the equation above the prime ’ denotes the derivative 
with respect to the conformal time and is the speed of 
sound of the fluid. For the present work H = 0 since it is 
a quantity that depends on the entropy perturbation and 
we will be focused with adiabatic perturbations only. In 
this regime the equation ([^ becomes homogeneous. In 
addition, the variables u and 9 are the Mukhanov-Sasaki 
variables defined through the following relations 

Po 

Po +Po 

By rewriting Eq. in such a way that rt is a function of 
the cosmic scale factor a and by considering a plane-wave 


4> 


AirGy/ Po + Po’ 



7 = 4'^- (9) 

with 0 < 7 < I for Eqs. The sound speed is = 
1 — 7 > 0 , and by use of we can write 

cl = 1-44. ( 10 ) 


Eor the following tachyon fields potential 


Viejr) 


Ay/l - 7o 

sinh2^<^’ 


( 11 ) 


some exact solutions of the field equation ([^ can be 
found, by assuming a linear dependence of the tachyon 
field with the cosmological time of the form 


4 — 4ot, 44 — 4o — 7o, 


( 12 ) 
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which is consistent with Eq. For this anzatz the 

Hubble parameter and its derivative are given by 


if=y|coth^^^t, 


H = - 


7oA 


sinh 


2 V37gA ^ 

' (13) 

which, after the integration gives the following cosmic 
scale factor 


a = ag 


sinh 



(14) 



It is customary to set Og = 1 for the actual scale fac¬ 
tor. The Eqs. (13) for the Hubble parameter and the 
derivative can be written in terms of the scale factor as 


_ 3 _ 1 

IP ~ “2^°TTa^ 



1 -t 

^370 


(15) 


The formulas just derived are valid for the standard 
tachyon. Nevertheless there exists two new kinds of ex¬ 
tended tachyon fields [21], [22|. In the last of those ref¬ 
erences it is shown that these two new types of tachyons 
could be derived from the standard tachyon field (0 < 
7 < 1) by use of a symmetry argument |26j . In the 
following we will consider the standard and the comple¬ 
mentary (1 < 7 ) tachyon field. 

The complementary tachyon field (pc is characterized 
by 1 < 0 c = 7 o and its potential can be deduced from 
the standard tachyon under the replacement 1 — 70 —t 
-(1 - 7 o) in Eq. ( 0 - resulting in 


V{Pc) 


AV 70 - 1 

sinh^ pc 


(16) 


By assuming again a linear dependence of the tachyon 
field with the cosmological time, Eq. (12), we can get 


the same Hubble parameter and its derivative, Eqs. (13) 


with 1 < 7 o. Note that, under these conditions and by 
use of Eq. (10), the sound speed becomes = 70 — 1. 
At first sight the energy density becomes 


V 


P(j> — 


Vl - 70 


and seems to have a divergence when 70 —1. Note 
however, that the potential 0 goes to zero in this limit, 
and the density remains finite. Thus the singularity 
7 o —>■ 1 is in fact an avoidable one. 

In Fig. the scale factor a is plotted as a function 
of t for the ACDM model and the standard and com¬ 
plementary tachyon field with different 70 values. The 
standard and complementary tachyon fields for values of 
7 o ~ 1 (70 = 0-99 and 70 = 1 . 01 ) are very close to the 
ACDM model. For these two values one infers from Fig. 
[ 1 ] that the cosmic scale factors increase with time, and 
there is no sensible difference between the solutions con¬ 
cerning the time evolution of the scale factors with the 


FIG. 1: Plot of cosmic scale factors as functions of time for the stan- 
dard (70 = 3/4 and 70 = 0.99) and the complementary (70 = 1.01 and 
7 o = 1.9) tachyon fields. 


concordance model. The two tachyon fields here behave 
as matter fields which are responsible for the decelerated 
regime. 

The standard tachyon field with 70 = 0.75 plays the 
role of an inflationary field which decays into matter. 
The analyses of Fig. shows that the evolution in time 
of the scale factor is less accentuated than that for the 
concordance model. 

On the other hand, the complementary tachyon field 
with 7 o = 1.91 that is close to represents stiff matter 
(70 = 2 ), the increase with time of a is much more ac¬ 
centuated than the one for the ACDM model. For a 
detailed analysis see m- 


IV. POWER SPECTRUM ANALYSIS 


We now proceed to present our results. Our aim is 
to compare the power spectrum results obtained for the 
tachyonization model for different values of 7 with the 
fiducial ACDM cosmology. 

In the following, we shall always refer to the dimen¬ 
sionless power spectrum A‘^{k,a) defined as 

a) = 16p{p + p) I U]^(a) p P, (17) 

with is the Fourier amplitude of the Mukhanov vari¬ 
able u, which is described by the Eq. (|^. It is therefore 
important to set up the initial conditions for the simula¬ 
tion. The following initial conditions for the variable 

m, m 


Wk(a) 


i 

V ^ A : 3/2 



a2i7fci/2’ 


(18) 


will be chosen. The reason is these are the minimal fluc¬ 
tuations allowed by the uncertainty relations. According 
to the current inflationary theory, the primordial quan¬ 
tum fluctuations are responsible for the large scale struc¬ 
ture formation. This election is the one which minimize 
the energy of these fluctuations. 
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We investigated the power spectrum for the tachy- 
onization model solving Eq. Q numerically with the 
initial conditions (18). These power spectrum was com¬ 
pared with the fiducial ACDM model corresponding to 
the following cosmological parameters: flmo = 0-30 for 
the total matter content, a cosmological constant con¬ 
tribution specified by flyvo = 1 ~ ^mo ^nd the Hubble 
constant i?o = 71kms~^Mpc~^. These parameters are 
in agreement with the ones obtained in WMAP9 [21] and 
the Planck mission l29l. 



FIG. 2: Plot of the dimensionless power spectrum as function of k 
at a = 1, for standard (70 = 3/4 and 70 = 0.99) and complementary 
( 7 q = 1.01 and 70 = 2 ) tachyon fields. 


The Figure [2] de scribes the a = 1 dimensionless power 
spectrum Eq. (|17|) for the tachyonization of the ACDM 
model, for the standard and the complementary tachyon 
field with 70 = 3/4, 70 = 0.99 and 70 = 1.01, 70 = 
2 respectively, and for the concordance model. These 
curves are the numerical solution of the Eq. @ with the 
corresponding initial conditions (18) for the above cases. 

We should remark that the model with 7 = 1 is the 
one which better imitates the expansion history of the 
universe corresponding to the fiducial model. However, 
their power spectrum differs considerably. The most ev¬ 
ident feature about this figure is the large discrepancy 
between the power spectrum for all the curves. As it 
can be seen from these results, there is a strong depen¬ 
dence of the power spectrum with the barotropic index 
7 o. For example, by changing from a barotropic index 
7 o = 3/4 to 7 o = 0.99 the dimensionless power spectrum 
varies in several orders of magnitude. Comparison with 
the ACDM model reveals that all the curves are strictly 
below the concordance model, with the exception of the 
one with 70 = 2. The last one is strictly above. Figure 
also shows that the lowest curve corresponds to the one 
with smallest barotropic index (70 = 3/4) and increases 
when 7 o gets larger. 

The Figure describes the power spectrum for the 
complementary tachyon held for 70 = 1.82, 70 = 1.91 and 
7 o = 2 and for the ACDM model. It is observed that the 
curve of the complementary tachyon held with 70 = 1.91 
coincides, has the same shape and amplitude, with the 
one of the concordance model. It can be seen from Figure 



FIG. 3: Plot of the dimensionless power spectrum as function of k at 
a = 1, for the complementary tachyon fields with 70 = 1.82, 70 = 1.91 
and 7 o = 2 . 


[^ that there is a fractional diherence of the power spec¬ 
tra of the complementary tachyon with 70 = 1.91 from 
that of the ACDM model is very small. 
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FIG. 4: Plot of the fractional difference of the power spectra of the 
complementary tachyon with 70 = 1.91 from that of the ACDM model. 


From Figure and m, it can be inferred that the 
models with 70 = 0.99 and 70 = 1.01 closely resemble the 
ACDM model. But when density perturbations are taken 
into account Figure shows that the two models with 
that barotropic index differ considerably. This result is 
physically relevant as it implies that the large structure 
formation that these models predict are far from obser¬ 
vations. The only special case, as we indicated above, are 
the ones with barotropic index very close to 70 = 1.91. 

Our interpretation about the discrepancies given above 
is the following. For the fiducial model the speed of sound 
Cg in is zero. The same happens with the tachyon 
model with 70 —>■ 1 , since the speed of sound is = 1 — 
7 o, so the equations look similar. However, this equation 
involves the factor (15) which differs for both models. 
When 7 o —)■ 1 the tachyon factor is 


H 

IP 




3 


1 


2 1 -b ■ 
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Instead for the ACMD model this factor is known to be 


H 

IP 


-3 


1 

1 + 


(19) 


where a has a value that depends of the density param¬ 
eters of the model. By comparing both factors it is seen 
that there is a coefficient 1 /2 of discrepancy between the 
two models. This simple difference can introduce radical 
change of behavior, even in linear differential equations. 
To give an example consider an hypergeometric equation 


z{l — z)u" + [a — (b + c + l)z]u' — bcu = 0 . 


This equation is clearly linear. One of its elementary so¬ 
lutions is the hypergeometric function F{a,b,c, z). This 
solution is convergent in the boundary |z| = 1 when [30] 

3fJ(6 -I- c — a) <0. 


Suppose that b and c have a fixed value, say 6 -|- c = 2. 
If for instance a = 3, then the condition above shows 
that F{a,b,c, z) is convergent in \z\ = 1. However, if we 
divide a with a factor 2, then a = 3/2 and the condition 
above is not satisfied. Thus the solution F{a, b, c, z) has 
a divergent behavior at the boundary \z\ = 1 . 

The example given above shows that a simple numer¬ 
ical factor in a linear equation can modify radically the 
behavior of the solutions. Unfortunately we do not know 
the exact solution of our equation Q, but our numer¬ 
ical results are suggesting that a subtle effect of this 
type is the one which gives the discrepancy between 
the power spectrum of the concordance model and the 
tachyon model with 70 —>■ 1. On the other hand, for 
7 o = 2 the factor (15) becomes the same as for the con¬ 
cordance model, but the speed of sound Cg is not zero 
and therefore the third term in Eq. Q do not vanish, 
thus the power spectrum is also different. Note that we 
did not even take into account the non trivial value of 
a in the factor (19), which also may generate discrep¬ 
ancies. The power spectrum is better approximated at 
some value between 70 —>■ 1 and 70 = 2 , which we found 
numerically to be 70 1.91, as shown in Fig. 


V. STATEFINDER DIAGNOSTIC 

The statefinder diagnostic is a diagnos introduced in 
|23j . which allows to distinguish between several dark 
energy models. This diagnostic is constructed starting 
with the scale factor a(t) and its derivatives up to the 
third order. The statefinder pair {s,r} [3T] [3^ [33] is 
defined by 


a H 

aiJ3 “ H3 , 

( 20 ) 

r — 1 

( 21 ) 

3(9- 1 / 2 )’ 


where q is the deceleration parameter, defined hy q = 
—a/aH'^. The statefinder pair is a geometrical diagnostic 
in the sense that it is constructed directly from the space- 
time metric. 


In the situations in which the space time is a flat FRW 
universe, the standard ACDM model corresponds to a 
fixed point {r,s\\cDM = {IjO}. This value is indepen¬ 
dent on the parameter of ACDM model and the redshift 
0 . Hence, we apply the statefinder diagnostic to the ta- 
chionization model for the different barotropic index 79 . 
By plotting the trajectories in the r — s phase diagram, 
the distance of the model from the concordance model 
can be probed, as discussed in [33] . 

By use of Eqs. (13), (14) and its derivates, the decel¬ 
eration parameter and Eq. (20) can be expressed as 




3 

2 



i) + 


1 

(z -k 1)3 T'o -k 1 


( 22 ) 


r 



(z -k 1 ) 3 T'o 
( z -k 1)3 t'o -k 1 ■ 


(23) 


The remaining parameter s can be obtained by substitut¬ 
ing equation (23) in equation ( |^ . We avoid to rewrite 
it here. 


Let us analyze the parametric plots of the Statefinder 
diagnostic. The Fig. shows the trajectories of {r, q} 
phase plane for four different values of the parameter 79 . 
This figure suggest that all the trajectories evolve to a 
fixed point in the future. The evolution for the trajec¬ 
tories with 79 = 0.99 and 79 = 1.01 (close to dust) is 
similar and overly far from the evolutions with 79 = 3/4 
and 7 o = 1.91. The phase plane for {s,q} is showed in 
Figure]^ Here also, the trajectories tend to a fixed point 
in the future, that coincides with the fixed point of the 
ACDM model (with {s,g} = {1,-1}). As expected, the 
evolution of the trajectories for the different barotropic 
index tend the one corresponding to a cosmological con¬ 
stant, so this model could explain the late time speed-up 
expansion of the universe. Finally, the {r, s} diagram 
can be seen in Fig. In this phase plane also, the 
trajectories evolve to a stable point in future which is 
corresponding to concordance model. 

In the Table |T] the currents values of the parameters 
q, r and s are represented for the different values of the 
barotropic index analysed in the work, 79 = 3/4, 79 = 
0.99, 7 o = 1.01 and 79 = 1.91. The present-day value 
of the deceleration parameter q{z = 0) S [—0.62; —0.56] 
as stated in the WMAP-9 report |3S], is below form the 
ones found in this work for all the barotropic index. In 
particular, for 79 = 1.91 the model is not in a acceleration 
phase, but goes to one in the future. The current values 
r(z = 0 ) and s(z = 0 ) in the model differ from 1 and 0 
respectively, for all the cases analyzed, see Table |T| 
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q 

FIG. 5: Trajectories in {t, g} phase plan with 70 = 3/4, 70 = 0.99, 
7o = 1.01 and 70 = 1.91. The black dot is the stable state of {r, q} in 
the future. 



FIG. 0! Trajectories in {s,q} phase plan with 70 = 3/4, 70 = 0.99, 
7 q = 1.01 and 70 = 1.91. The black dot is the stable state of {s, q} in 
the future. 


70=1-91 ; 

\ 


7o=1.01 

k ACDM 

70=0.99 

\ 

\ 70=3/4 


FIG. 7'. Trajectories in {r, s} phase plan with 70 = 3/4, 70 = 0.99, 
7o = 1.01 and 70 = 1.91. The black dot is the stable state of {r, s} in 
the future. 


Tachyonization of the ACDM model 

70 

q{z = 0) 

r(z = 0) 

s(z = 0) 

70 = 3/4 

-0.44 

0.58 

0.15 

7o = 0.99 

-0.26 

0.98 

0.01 

70 = 1.01 

-0.24 

1.02 

-0.01 

70 = 1.91 

0.43 

4.91 

-19.31 


TABLE I: The present values of the currents q, r and s 
parameters for different barotropic index. 


VI. SUMMARIES AND CONCLUSION 


In the present work the power spectrum for the tachyon 
models of m. m was investigated by numerical meth¬ 
ods, by taking (18) as initial conditions. The power 
spectrum of the standard and complementary tachyon 
was compared with the spectrum of the fiducial ACDM 
model, since we know that the standard ACDM model 
is pretty much consistent to the current cosmological ob¬ 
servations. This was done by taking into account cos¬ 
mological parameters which are in agreement with the 
ones obtained in WMAP9 [3S] and the Planck mission 
[29] . It was found that both spectrums are quantitatively 
different and can be distinguished, even for small A val¬ 
ues. The power spectrum obtained strongly depends on 
the value of the barotropic index 70 of the tachyon in 
consideration, although the qualitative behaviour of the 
power spectrum curve is similar for all the cases. For in¬ 
stance, the power spectrum for the cases with 70 = 3/4 
differ in several orders of magnitude from the one with 
7 o = 0.99. We have found that the tachyon model that 
better agrees with the ACMD is the one with barotropic 
index 70 = 1,91. This corresponds to a complementary 
tachyon. The curves for the power spectrum correspond¬ 
ing to 7 o < 1,91 are strictly below the ACMD model, 
while the opposite holds for 1.91 < 70 < 2. 


In addition, the statefinder diagnostic for the minimal 
case has also been studied. The plot of the phase plan 
of the pairs {r,q}, {s,^} and {r, s} suggested that the 
trajectories of the all cases, for the first pair, tend to the 
state which has a non vanishing distance from ACDM 
model. In contrast, the last two pairs reaches a stable 
state which is corresponding to the concordance scenario, 
so the tachyonization of the ACDM model is capabie to 
explain the late time accelerating phase of the universe 
expansion. In the Table [J it is shown the deceleration 
parameter q for different values of the barotropic index, 
and it is found that for 70 = 1.91 the g(0) = 0.43, which 
implies that the universe is decelerating instead of ac¬ 
celerating. An odd observation is that this value of 70 
gives a similar power spectrum that the fiducial model, 
although a decelerating universe is ruled out by the ob¬ 
servations. On the other hand, the values of 70 ~ 1 which 
can result to similar expansion history as ACDM model, 
as the Figure [ 1 ] shows, have several orders of magnitude 
difference in amplitude of the power spectrum in compar¬ 
ison to the fiducial model. So it seems that the assumed 
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Tachyonic model cannot explain all observations simul- proved the description of the results. The authors are 
taneously. supported by the CONICET, Argentina. 

Our work did not considered inflationary aspects of 
this tachyon fields, and it may be of interest to study 
them. We leave this for a future investigation. 
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